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P_( ' Abstract 

' The relaxation limit in critical Besov spaces for the multidimensional compress- 

• ible Euler equations is considered. As the first step of this justification, the uniform 

, (global) classical solutions to the Cauchy problem with initial data close to an equi- 

librium state are constructed in the Chemin-Lerner's spaces with critical regularity. 
Furthermore, it is shown that the density converges towards the solution to the 
porous medium equation, as the relaxation time tends to zero. Several important 
' estimates are achieved, including a crucial estimate of commutator. 

^ : 

I Keywords: compressible Euler equations, classical solutions, relaxation limit, Chemin-Lerner's 

ff^ ■ spaces 

o\ 

O ; AMS subject classification: 35L25, 35L45, 76N15 



1 Introduction and Main Results 

In a suitable nondimen 
polytropic fluid read as 



X 

' In a suitable nondimensional form, the multidimensional compressible Euler equations for a 



r dtp + V ■{p^r) = 0, 

\ dtipv) + V • (pv v) + VP = ^ > 

Here p = p{t, x) is the fluid density function of {t, x) € [0, +cxd) x with d > 2; ^ = {v^ ,v'^ , ■ ■ 
■,v'^)^{T represents the transpose) denotes the fluid velocity. The pressure P = P{p) satisfies 
the usual 7-law: 

P{p) = Ap'<{^ > 1) 



* E-mail : j iangxu_79 @yahoo . com. cn 
^^E-mail: wangzejun@gmail.com 



1 



where ^4 > is some physical constant, the adiabatic exponent 7 > 1 corresponds to the 
isentropic flow and 7 = 1 corresponds to the isothermal flow; < r < 1 is a (small) relaxation 
time. The notation V,® are the gradient operator and the symbol for the tensor products of 
two vectors, respectively. 

In this paper, we are going to study the Cauchy problem of the compressible Euler equations 
(II. ip subject to the initial data 

(p,v)(0,x) = (po,vo). (1.2) 

Our first interest is, for fixed r > 0, to investigate the relaxation effect on the regularity and 
large-time behavior of classical solutions. As shown by \19\ 121] . if the initial data are small in 
some Sobolev space H^{M.'^) with s > 1 + d/2 (s G Z), the relaxation term which plays the role 
of damping, can prevent the development of shock waves in finite time and the Cauchy problem 
(|l.ip - (|1.2p admits a unique global classical solution. Furthermore, it is proved that the solution 
in [19] has the L°° convergence rate (l + t)~^/^((i = 3) to the constant background state and the 
optimal L'P{1 < p < 00) convergence rate (1 + t)~'^/2(i-i/p) in general several dimensions [21], 
respectively. For the one-dimensional Euler equations with relaxation, the global existence of a 
smooth solution with small data was proved by Nishida |17j . and the asymptotic behavior of the 
smooth solution was studied in many papers, see e.g. the excellent survey paper by Dafermos 
[3] and the book by Hsiao [8j. In addition, for the large-time behavior of solutions with vacuum, 
see PllO]. 

Another main interest is to justify the singular limit as r — )• in (jl.ip . First, we look at the 
formal process. To do this, we change the time variable by considering an "0(l/r)" time scale: 

(/,^vO(.,x) = (p,v)(^,x). (1.3) 

Then the new variables satisfy the following equations: 



(1.4) 



with initial data 



(p^v-)(x,0) = (po,vo). (1.5) 



At the formal level, if we can show that ^ ^ is uniformly bounded, it is not difficult to see that 
the limit of /o"^ as r — >• satisfies the porous medium equation 

dsM - APiM) = 0, 

which is a parabolic equation since -P(AA) is strictly increasing. 

This singular limit problems for hyperbolic relaxation to parabolic equations have attracted 
much attention. By means of compensated compactness theory, Marcati and his collaborators 
[I3l[l4l[l6] systematically studied this diffusive limit of generally quasi-linear hyperbolic system, 
also including the present Euler equations (jl.ip for weak solutions. When 7 = 1, Junca and 
Rascle verified the convergence of solutions to the isothermal equations (II. ip towards the 
solution to the heat equation for arbitrarily large initial data in BV(M) that are bounded away 
from the vacuum. Coulombel and Goudon [3] fell back on the classical energy approach and 
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constructed the uniform smooth solutions to the isothermal Euler equations and justified the 
relaxation limit in some Sobolev space H^(M.'^){s > 1 + d/2, s G Z) (in x). 

In the present paper, we will improve Coulombel and Goudon's work [3J such that the 
result may hold in the critical space with the regularity index a = 1 + d/2 (a larger space). 
Indeed, we choose the critical Besov space B2i{M.'^) in space- variable x rather than H'^(M.'^) 
as the functional setting, since i?2i(I^'^) is a subalgebra of W^'°°. Starting from this simple 
consideration, based on the Littlewood-Paley decomposition theory and Bony's para-product 
formula, we first construct the (uniform) global existence of classical solutions in the framework 
of the Chemin-Lerner's spaces L^(i?p ,.) in [2j, which is a refinement of the usual spaces L^{Bp j.). 
Then, using Aubin-Lions compactness lemma, we perform the relaxation limit of (jl.ip - (jl.2p in 
Besov spaces. 

Main results are stated as follows. 

Theorem 1.1. Let p > be a constant reference density. Suppose that po — p and vq S 
i?2 i(IK'^)(<7 = 1 + d/2), there exists a positive constant Sq independent of t such that if 

IKpo - ^,vo)||Ba^(Rd) < 6o, 

then the Cauchy problem U.1\ }- [TI^) has a unique global solution (p, v) satisfying 

(p,v) G C\R+ X M"^) 

and 

[p -p,^)e C(i?J,i(M'^)) n C\B^\R'')). 

Furthermore, the uniform energy inequality holds 

Up- 

+Ao{ 



1 



r 



+ 



V^Vp 



< CoWipo 



P, ^o)\\bi^( 



(1.7) 



where < r < 1, Aq and Cq are some uniform positive constants independent of t. 



Remark 1.1. In comparison with that in [3j, Theorem ll.il depends on the low- and high-frequency 
decomposition methods rather than the classical energy approach. As shown by ourselves 
the low-frequency estimate of density for the Euler equations (jl.ip is absent. Then, we overcame 
the difficulty by using Gagliardo-Nirenberg-Sobolev inequality (see, e.g., 0) to obtain a global 
classical solution, however, the result fails to hold in the critical Besov spaces mentioned above. 
To obtain the desired result, in the current paper, we add the new context in the proof of global 
existence. Indeed, some frequency-localization estimates in Chemin-Lerner's spaces are devel- 
oped, including a crucial estimate of commutator, for details, see Proposition 14.11 Proposition 
16.11 and Corollary 16.21 

Based on Theorem ll.lt using the standard weak convergence method and Aubin-Lions com- 
pactness lemma in |18], we further obtain the relaxation limit of ()l.ip - ()1.2p in the larger frame- 
work of Besov spaces. 
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Theorem 1.2. Let (/O, v) be the global solution of Theorem \l.l\ Then 

— p is uniformly bounded in C(M'*', i?2 

is uniformly bounded in L^(M+, ^2 

r ' 

Further, there exists some function TV G C(M+,n + B^^{W^)) which is a global weak solution of 

For any < T,R < oo, {p'^{s,x)} strongly converges to N{s,x) in C([0, T], (S2l^(-Sr-)) 

as r — 7- 0, where 5 G (0, 1) and denotes the ball of radius r in M*^. In addition, it holds that 

||(A/'(S,-) -pllB-^(Md) < C'oIKPO -P,Vo)||iJ-^(l;d), s > 0, (1.8) 

where C'q > is a uniform constant independent of t. 

Remark 1.2. Compared with that in [3], the relaxation convergence of classical solutions holds in 
the Besov spaces with relatively lower regularity. To the best of our knowledge, this is the first 
result for the Euler equations (jl.ip in this direction. Therefore, Theorem 11.21 gives a rigorous 
description that the porous medium equation is usually regarded as an appropriate model for 
compressible inviscid fluids. In addition, let us also mention that the limit result is generalized 
to be true for general adiabatic exponent 7 > 1 but not the only case 7 = 1 in [3]. 

The paper is organized as follows. In Section 2, we briefly review the Littlewood-Paley 
decomposition theory and the characterization of Besov spaces and Chemin-Lerner's spaces. In 
Section 3, we reformulate the equations (jl.ip as a symmetric hyperbolic form in order to obtain 
the effective frequency-localization estimate and present a local existence result for classical 
solutions. In Section 4, using the high- and low-frequency decomposition methods, we deduce 
the frequency-localization estimate in Chemin-Lerner's spaces, which is used to achieve the 
global existence of uniform classical solutions. Section 5 is devoted to justify the relaxation 
limit for the Euler equations (jl.ip . Finally, the paper ends with an appendix, where we give the 
proof of estimates of commutator. 

Notations. Throughout this paper, C > is a generic constant independent of r. Denote 
by C([0,T],X) (resp., C^{[{),T],X)) the space of continuous (resp., continuously differentiable) 
functions on [0, T] with values in a Banach space X. We often label ||(a, 6, c, (i)||x = + 
ll^llx + ||c||x + Mllx) where a,b,c,d £ X. Here and below, we omit the space dependence for 
simplicity, since all functional spaces are considered in W^. Moreover, the integral f^^ fdx is 
labeled as J / without any ambiguity. 

2 Preliminary 

For convenience of reader, we try to make the context self-contained, in this section, we briefly 
review the Littlewood-Paley decomposition theory and some properties of Besov spaces and 
Chemin-Lerner's spaces. For more details, the reader is referred to p4i5j. 

Let {(f, x) be a couple of smooth functions valued in [0, 1] such that ip is supported in the 
shell C(0, f , I) = G M'^ll < ICI < f }, X is supported in the ball B(0, |) = G M'^H^I < 5} 
and 

00 

g=0 
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Let S' be the dual space of the Schwartz class S. For / G S', the nonhomogeneous dyadic blocks 
are defined as follows: 

A_i/ := x{D)f = oj*f with oj = 
A J := ip{2-'iD)f = 2^^^ J uj{2'^y)f{x - y)dy with a; = J"" V, if > 0, 

where * the convolution operator and J-~^ the inverse Fourier transform. The nonhomogeneous 
Littlewood-Paley decomposition is 

/ = 5^ A,/ V/ e S'. 

q>-l 

Define the low frequency cut-off by 

S,f := A,/. 

p<q-l 

Of course, Sof = A_i/. Moreover, the above Littlewood-Paley decomposition is almost orthog- 
onal in L^. 

Proposition 2.1. For any f G 5'(IR'^) and g £ S'(M.'^), the following properties hold: 

ApAgf = if \p-q\>2, 

Ag{Sp-ifApg) = Q if \p-q\>5. 

Having defined the linear operators Aq{q > —1), we give the definition of Besov spaces and 
Bony's decomposition. 

Definition 2.1. Let 1 < p < oo and s G M. For 1 < r < oo, Besov spaces Bp ,,. C S' are defined 
by 

f e Bl, ^ =: ( ^ (2'?^||A,/|U.f ) ' < oo 

q>-l 

and Bp^^ C 5' are defined by 

f e S^,oo ^ II/IIbIoo =■■ sup 2''^||A,/|U. < oo. 

q>-l 

Definition 2.2. Let f,g be two temperate distributions. The product / • g has the Bony's 
decomposition: 

f-g = Tfg + Tgf + R{f,g), 
where Tfg is paraproduct of g by f, 

Tfg= Ap/A,5 = ^Vi/A,i; 

p<q-2 q 

and the remainder R{f, g) is denoted by 
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As regards the remainder of paraproduct, we have the following result. 
Proposition 2.2. Let (si,S2)€M^ and 1 < p,pi,p2,r,ri,r2 < oo. Assume that 
111 111 

-< \ <1, -< \ , and si + S2>0. 

p Pi P2 r ri r2 



Si+S2+rf(- -) 

Then the remainder R maps -Bp J x Bp^ ^.^ in Bp^r ^ '"^ CLud there exists a constant C 
such that ^ ^ 

< -^^^\\f\\B;ij9\\B;i.,- 

Some conclusions will be used in subsequent analysis. The first one is the classical Bernstein's 
inequality. 

Lemma 2.1. Let A; G N and < Ri < i?2- There exists a constant C , depending only on Ri, R2 
and d, such that for all 1 < a <b < 00 and f £ L^, 

Supp Tf C B(0,i?iA) sup Wd'^fh^ < C'=+iA'=+^(^-^)||/||ia; 

\a\=k 

Supp ^/ C C{0,RiX,R2X) C-^-'xHfh^ < sup Wd^fh^ < c^+^x'^WfU'^. 

\a\=k 

Here J- f represents the Fourier transform on f . 

As a direct corollary of the above inequality, we have 
Remark 2.1. For all multi-index a, it holds that 

||5"/b.„<C||/||^.+H. 

The second one is the embedding properties in Besov spaces. 

Lemma 2.2. Let s G M and I < p,r < 00, then 

B^„ ■-^ B^f whenever s < s or s = s and r < r: 

s-(i(i-4) 

Bp J. ^ Bp^ ^ whenever p > p; 

B^Jjil <p<oc)^Co, C n L~, 

where Cq is the space of continuous bounded functions which decay at infinity. 

The third one is the compactness result for Besov spaces. 

Proposition 2.3. Let l<p, r<oo, sGM and e > 0. For all 4> G C^iW^), the map f ^ (j)f 
is compact from B^p{R'^) to B^ ,.{R'^). 

On the other hand, we also present the definition of Chemin-Lerner's spaces first incited by 

.s ^ 
p,rJ' 



J.-Y. Chemin and N. Lerner [2], which is the refinement of the spaces L^{B!^ 
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Definition 2.3. For T > 0,s < r,9 < oo, set (with the usual convention if r = oo) 

9>-l 

Then we define the space Lj^{Bp j,) as the completion of S over (0, T) X by the above norm. 
Furthermore, we define 

Ct{b;^,.) ■.= L^{B;^,.)nc{[o,T],B;^,.) 

and 

:= {/ G C\[0,T],B;^,)\dtf G 
The index T will be omitted when T = +00. Let us emphasize that 
Remark 2.2. According to Minkowski's inequality, it holds that 

Then, we state the property of continuity for product in Chemin-Lerner's spaces L^{Bp j.). 
Proposition 2.4. The following estimate holds: 

whenever s > 0,1 < p < 00, 1 < 9,9i, 02,0-^, 64^ < 00 and 

1 _ 1 1 _ 1 1 
6 61 02 O3 Oi 

As a direct corollary, it holds that 

||/5llz«(B.,j<C^II/llz^i(B|,.)ll^^llz^^(B|,J 
whenever s > d/p, ^ = + ^■ 

In addition, the estimate of commutators in Lj,{Bp ^) spaces is also frequently used in the 
subsequent analysis. The indices s,p behave just as in the stationary case [U [5] whereas the 
time exponent p behaves according to Holder inequality. 

Lemma 2.3. Let 1 < p < 00 and 1 < 9 < 00, then the following inequality is true: 
2^^||[/,A,M5ll4(Lp) <Cc,||/||~.,^^.^^||<7||£..(^.^), s = l + d/p, 

where the commutator [•, •] is defined by [f,g] = fg — gf , the operator A = div or V , C is a 
generic constant, and Cq denotes a sequence such that ||(cg)||;i < 1, | = ^ + ^• 

Finally, we state a continuity result for compositions to end up this section. 

Proposition 2.5. Let s > 0, 1 < p,r, 9 < 00, F e M) with F{0) = 0, T e (0,oo] 

and V £ L^TiB'p r) n L^{L°°). Then 



7 



3 Symmetrization and local existence 



In terms of the ideas in |19j . we introduce a new variable (sound speed) which transforms the 
equations (jl.ip into a symmetric hyperbohc system. For the isentropic case (7 > 1), denote the 
sound speed by 

and set = ip{p) corresponding to the sound speed at a background density p > 0. Let 

2 



-y(^(p)-^)- 



7- 

Then the Euler equations (jl.ip is transformed into the symmetric form for classical solutions: 



dtg + V'divv = — v • Vg — ^^^gdiw, 
dtv + i^Vg + iv = -V • Vv - ^gVg. 



2 ^;iV' (3.1) 

(v^i/yvy^-v — — V - vV 2~g^ g- 

The initial data (jl.2p become 



ig,v)\t=Q = {gQ,vo) (3.2) 

with 

2 

00 = 7(^(Po) - V')- 

7-1 

Remark 3.1. The variable change is from the open set {(/j, v) G (0, +00) x M'^} to the whole 
space {(£>,v) G RxR^}. It is easy to show that for classical solutions (p, v) away from vacuum, 
(fLU-dllll) is equivalent to (I3T])-(I32|) with ^g + > 0. 

Remark 3.2. For the isothermal case 7 = 1, let us introduce the enthalpy change g{t,x) = 
^/A{lnp — In^o). In this case, the equations (jl.ip can be transformed into the system ()3.ip with 
7 = 1, the reader is referred to [7] for more details. In subsequent sections, we focus mainly on 
the case 7 > 1, since the isothermal case can be dealt with at a similar manner. 

Recently, we have achieved a local existence theory of classical solutions in the framework of 
Chemin-Lerner's spaces for compressible Euler-Maxwell equations, see [22j. Actually, the new 
result is applicable to generally symmetrizable hyperbolic systems, including the current Euler 
equations of special form. Here, we present the result only and the details of the proof are 
omitted for brevity. 

Proposition 3.1. For any fixed relaxation time t > 0, assume that (£>0)Vo) G B21 satisfying 
■^^£"0 + ip > 0, then there exists a time Tq > (depending only on the initial data) and a 
unique solution {g,v) to [3l]) - [3^) such that {g,v) G C^([0, Tq] x R'^) with ^g + ip > for all 
t G [0,To] and (^>,v) G CroiBl^) DCUB^T^). 



4 Global existence 

In this section, we first establish a crucial a priori estimate in Chemin-Lerner's spaces. Then 
by the standard boot-strap argument, we obtain the global existence of classical solutions of 

The a priori estimate is comprised in the following proposition. 
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Proposition 4.1. Let (£», v) e Ct{B^ i) nC^{B^^'^) be the solution of ^31^-^3^ for any given 
time T > 0. There exist some positive constants 5i,Xi and Ci independent of t such that if 



(4.1) 



then 



+Ai 



+ 



< Ci\\{eo,^ro)\\Bl,. 



(4.2) 



Proof. The proof of Proposition 14.11 in fact, is to capture the dissipation rates of (£», v) in turn 
by using the low- and high-frequency decomposition methods, so we divide it into several steps. 

Step 1. The Lf{B^^^) estimate of (^,v) and the L\{B^^^) one of v; 
Firstly, we complete the proof of step 1. Applying the localization operator Ag to p.ip 
yields 



dt\Q + ^A^divv + (v • V)Ag^ = [v, Ag] • - \iAg(^divv), 



9t AgV + ijj/^qSjQ + (V • V) AgV + 



[v,A,].Vv-3:2^A,(£.V£>), 



(4.3) 



where the commutator[-, •] is defined by [/,(/] = fg — gf- 

Multiplying the first equation of (|4.3p by AgQ and the second one by A^v respectively, then 
integrating them over M*^, we get 



ld_ 
2'dt 



|A,^,||^2 + ||Agv||^2)+i||A,v||i2 



i j divv(|A,^|2 + |Aqv|2) + y"{[v, AJ • Vg/^qQ + [v, Ag] • VvA^v} 

+ / ^qQi^Q ■ ^qV) + / [q, Aq]VQ ■ AgV + / [q, Ag]divvAg^. (4.4) 



In what follows, we first bound the low frequency part of the quality (14. 4p . By performing 
integration by parts and using Holder- and Gagliardo-Nirenberg-Sobolev inequalities, we have 
{d > 3) 



i|(||A_,,||i. + ||A_,v 



lL2 



+ -||A_iv 



< ||v||j^d||A_i^|| 2d ||A_i V<)||j;2 + ||Vv|| TOO ||A_i v||^:3 

+ ||[v,A_i] •V£*||^^||A_i£.||^^ + ||[v,A_i]- Vv||i2||A_iv||i2 



7-1 

2d_ II A_iv||i2 H —II [q, A_i]V£»||i2 II A_iv||i2 



+^l|V.k.||l|A-,.„^^ 

+1^11 [^,, A_i]divv||^^ II A_i£.||^^ 
< ||v||^d||A_iV£)||22 + ||Vv||l^||A_iv||22 + ||[v, A_i] • V^|| 2d ||A„iV£>||i2 
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+ ||[v,A_i] -Vvll^allA^ivllia + 



7-1 



||V^>||i.||||A_iV^>||i2||A_iv||i2 



— ||[£),A_i]divv||^^||A_iVi?b2. (4.5) 



Integrating ()4.5p with respect to t G [0, T] implies 



< 



- ||A_i£,||i. + ||A_HliO 



r 



l|v||L2(^d)||A_lV£.||i2(i2)||A_lV£.||ioo(i2) + ||Vv||ioo(ioo)||A_lv||22^^2^ 

+ ||[v, A_i] ■ Vi?||^^^^^^||A_iV£>||i2(^2) + ||[v, A_i] ■ Vv||^2(i2)||A_iv||^2(^2) 



7-1 



7-1 , 
2 ' 
7-1, 



l|V^'||LJ>o(id)||A_lV£»||i2(i2)||A_lv||^2(i2) 
A_l]V£<||i2(i2)||A_lv||i2(i2) 
[£>,A_l]divv|| 2d ||A_lV^||i2(i2). 



2 "Li(L^y 

Then multiplying the factor 2~'^^ on both sides of (j4.6p . we can get 



(4.6) 



< 



1 



|A_i£>||i2 + ||A_iv||i2 + ||A_iv 



-2(T 



-2o- 



^iqqWI^ + ||A_ivo||i2 +Cc 



liv^ll 



L5?(BJ-1) 



|Vv| 



^oo^R<T-l^||v||~2 ^ + C'c_;^||v|L 



|Vv||r2,„._i 



L2(B27i)ll'^llL-(B2-i) 
l^^?IIZ-(B-7l)ll^^'IIZ2.o.-l^l|V||r 



l^?llL-(BJjll^^llL?.{iJ--i)l|v|lL2{B,-,,) 



|Vf?|| 



-'^T (-^2,1 ^) -^tC-S?,!) 



l|V£>||r: 



(4.7) 



where we used Remark 12.21 Lemma [2.31 and Corollary 16.21 which will be shown in the Appendix. 
Here and below C > denotes a uniform constant independent of r; {c_i} denotes some 
sequence which satisfies ||(c_i)||;i < 1 although each {c_i} is possibly different in (|4.7|) . 

Next, we turn to estimate the high-frequency part (g > 0) of the quality (|4.4p . With the aid 
of Cauchy-Schwartz inequality, we have 



^|(||A,,||i. + ||A,v||i.)+i||A,v||i. 



2dt 

< ^\\Vv\\L^{\\A,e\\h + IIA,v||io + ||[v. A,] • VeWvAWQh^ 



+ ||[v,A,]-Vv||i2||A,v||i2 + 



7-1 



||V£>||l^||A,£.||^2||A,v| 



L2 



+l^\\[g,Ag]Vgh2\\Agvh2 + ^ 



||[£),A,]divv||i2||A,^||i2. 



(4.8) 
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By integrating (j4.8p with respect to t G [0,T], we arrive at 



^\\\^<iQ\rL2 + IIA.vlli^j ^ + -||A,v||^2(^2) 

+ 11 [V, A,] • V£)||i2(i2) II A,£>||i2(i2) + II [V, A,] • VV||^2(^2) II AgV||^2(^2) 

7 - 1 

+ ^— l|V£'||L-(L-)l|Ag£'||i2(i2)||AgV||^2(i2) 

II [^', Ag]divv||^2(^2) II A,^^||i2(^2). (4.9) 
Then multiplying the factor 2^'^°^ on both sides of (14. 9p and using Lemma 12.3^ we obtain 

-22«'^(||A,£,||i. + ||A,v||i.) + --||A,v||2,(^,^ 
< i22«'^(||A,^,o|li2 + l|A,vo||i2) +Cc2||Vv||£.(^.-.)(||V^,||j^(^.-i)||^>y^ 

+Cc^||v||j^(^^,^^^||v|||2^^^^,^^^ + Cc^llV^II;^,,^^^^^^ 
+'^c2||V^>||j^(^.-l)||V£>||j2^(^.-l)||v||j2(g^ 

+'^c2||Ve||;^^(^.-i)||v||j2(^^.^^)||Vf5||j2(^.^ (4.10) 



< ^l|Vv||i2(^^)(||Ag^||^2(^2)||Ag£.||iC«(^2) + ||AgV||^2(^2)||AgV||iC«(i2)) 



where we have used the fact ||AgV/||^2 25||Aq/||2,2(g > 0) derived by Lemma [2.11 The 
constant C > is a uniform constant independent of r; {cq} denotes some sequence which 
satisfies ||(cq)||;i < 1 although each {cg} is possibly different in ()4.10p . 
To conclude, combining (|4.7p with (|4.10p gives 



1 , . nzqa 

-22'^'^(||A,£.||i. + ||A,v||i.) + --||A,v||^.(^,) 
< i22''-(||A,£.o|li2 + l|A,vo||i2) +Cc^||v||j.(^.^)(||V^>||£^(^.-.)||£,||j^ 

+<^''lyh^{BlJ\'^(^hl(B-,-')M^^^ (g > -1). (4.11) 

By Young's inequality, we have 

2^"(l|Ag£^||L5?(L2) + ||AgV||i^(^2)j + ^fi^||AgV||^2(i2) 



< 



C2'''^(||A,£,oIIl^ + l|A,vo||i2 j + Cc,^||^||j^(^.^) (^_||v||~, ^^^^^^^ + V^\\Ve\\i2^^s.-i) 



+Ccg^ML¥(Bi,)-^MLUBi,) il > -1)' (4-12) 
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where /ii is a positive constant independent of r. 

Summing up (14.120 on q > —1, we immediately get 



< 



1 



+C^I|v||Zoo(^.^)-^||v||j;.(^.^). (4.13) 

step 2. The L^B^^) estimate of V^i. 

Using the second equation of ()3.ip . we have 



^/;V£»= -(^atV+ -V + V- VV+ ^^y^^V^I^. (4.14) 

Apply the operator Ag to (I4.14p to get 

^TAgdtV + AqV - r[v, Aq]Vv + TV • AgVv 



^ ^T[Q,Ag]Vg+^^TgAgVQ). (4.15) 



Integrating the resulting equality over R*^ after multiplying AqVg, we have 

i;T\\AgVQ\\l2 

^rAgdtv + Agv - t[v, Ag]Vv + rv • A^Vv 



2 

where the first integral can be estimated as 



^ ^T[g, Ag]V0 + I-ItqA^Vq) ■ AgVg, (4.16) 



- T y Agdtv ■ AqVg = y AqdivajvAqgi 

= J AqdivvAqQ -T j AqdivvAqdtg 
= J AqdivvAqQ 

—T I ArtdivvAof — ^idivv — v • — — ^>divv 



< r— / AgdivvAg^ + T^lldivvll^a +r||v||L^||Agdivv||i2||AgV^?|U2 
+T\\Aqdiy^L2\\[v,Aq]Vg\\L2 + ^T\\0U^\\d\vv\\l, 
+ ^T||divv||i2||[^?,A,]divv||i2. (4.17) 



12 



Remark 4.1. In the inequality (j4.17p . the information behind the mass and momentum equations 
of ()3.ip help us eventually to estimate the term 9fV well. Otherwise, as in we have to establish 
an auxiliary inequality with respect to the variable {Qt,'Vt) to close the a priori estimate, which 
leads to the tedious proof of global existence consequently. 

Together with (j4.16p - (|4.17p . we are led to the estimate 
d f 

< ^J^J AgdivvAg£) + TV5||divv||i2 + ||Agv||i2||A,V£)||i2 
+r||v||L^||A,divv||i2||A,V^||i2 

+r||Agdivv|U2||[v, A,]V^?|U2 + ^^t\\q\\l^ \\Agdivv\\l2 



+^^t|| A,divv||i2 II [g, Ag]divv||i2 + t||v||l^ || AgVv||i2 1| A,V^>||i2 
+r|| [v, A,]Vv||^2 WA.Veh^ + ^II^^IIl^ II A,V^||i. 



+^T||[i,,A,]V^||i2||A,V^>||,.2. (4.18) 



Integrating (|4.18p in t G [0, T] gives 

< r(||A,divv||i2||Ag^||i2 + ||Agdivvo||L2||A,^o||L2) +rV;||A,divv||22^^2^ 

+ l|Agv||i2(i2)||A,V£>||i2(i2) +r||v||io=(i^)||A,divv||i2^(i2)||A,V£.||i2^(i2) 

7 — 1 

+T||A,divv||i2^(i2)||[v, A5]V£>||i2(i2) + T||£i||iop(ioo)||Agdivv||^2^^2^ 
7-1 

+^— ^l|Agdivv||i|(i2)||[^, Ag]divv||^2(^2) +T||v||ioc(ioo)||AgVv||i2^(i2)||AgV£»||i2(^2) 

+t||[v, Ag]Vv||i2^(^2)||A,V^>||i^(i2) + ^||^>||i5?(L-)||A,V^?||22(^2) 

+ ^T||[£,,A,]Vf?||^2(^2)||A,Vf?||^2(^2). (4.19) 

Multiply the factor 22'?('^-i) on both sides of (lil^ to get 

^22.(-i)||A,V,||i2(^2) 

< C'^c2(||divv||;^^(^.-i)||£>||j^(^.-l) + ||divVo||^.-l||f?o||B^._-l) 
+Crc^||divv|||^^^^._i^ + Cc2||v||;^2^(^.-i)||V£>||22^(^.-i) 
+CrCpv||2;^(^.-i)||divv||2;2(^.-l)||Vf?||£2(5^^^^ 

+Crcg||^»||j^(^^,^^)||v|||2 + Crc^||v||^^ 

+^^c'II^?IIl??(bj,)I|V^>II|^(s.-1), (4.20) 
where we used Lemma [231 and {cg} denotes some sequence which satisfies ||(cg)||;i < 1. 
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Then it follows from Young's inequality that 



Cc^. 



1 

1 



+^C,J||£.||j;^(^.^)V^||Vf?||j-,(^.-.), (4.21) 



where we have used the smallness of r(0 < r < 1). 

Finally, summing up (|4.2ip on (7 > —1, we deduce that 



1) 



, ) 

1 



+ (S2,i)77"''"4(S?,i) + ^\jML^{Bl,)^\\^e\\^(Blj^y (4.22) 



Step 3. Combining the above analysis. 

Combining with (I4.13P and (I4.22p . we end up with 



1 



+C'^l|v|lz-(Bjj^l|v||£2^(Bj^^) +C'if(||(£',v)||j^(^^.^^) + \\{qo,^o)\\bi, 

1 



+'^^^Ml^(bi,)^MlI{bi,) + CK^\\Q\\i^^j,.^^^^V^\\VQ\\i2^^^.-iy (4.23) 

where -fC > is a uniform constant independent of r . In order to eliminate the term 1 1 ( v) 1 1 j-o^ ^.^^ ^ 
and the singular one ||v||2;2 (^<t -^I Vt, we take the constant K such that 

Furthermore, it is not difficult to obtain 
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1 



1 



1 



1 



1 



< C\\{e,,^,)\\B^^^^+C^5r[-^\M^^^^^^^^ (4.24) 

where we have used the a priori assumption (j4.ip in the last step of (j4.24p . 
Lastly, we choose the positive constant 5i satisfying 

Cv^<min{^,K}, 

then the desired inequality (|4.2p follows immediately. □ 

With the help of the standard boot-strap argument, for instance, see [l5], Theorem [LT] fol- 
lows from the local existence result (Proposition 13. ip and a priori estimate (Proposition 14. 1^ . 
Here, we give the outline of the proof. 

Proof of Theorem If the initial data satisfy ||(£>0) vo)||_bj^ < by Proposition 13. H then 
we determine a time Ti > 0(ri < Tq) such that the local solutions of (j3.ip - (|3.2p exists in 
CxiiBni) and ||(i?, v)||roo s < 5i. Therefore from Proposition 14.11 the solutions satisfy the 

I ^Ti'^ 2,l) 

a priori estimate IK^*, v)|| jcx, ^ < Ci||(£»o, vo)||bj^ < ^ provided \\{qo,vo)\\bi^ < Thus 

by Proposition 13.11 the system (j3.ip - (|3.2p for t >Ti with the initial data {g,v)(Ti) has again a 
unique solution (£>, v) satisfying IKe*, v)||roc ^r<t n < Si, further ||(£>, v)||roo ^ < Si. Then 

by Proposition 14.11 we have IK^j, v)||7oo /r<t < C'iIK^jQ) vo)||_b;^, < ^. Thus we can continuous 

^2Ti^ 2,l) ^'^ ^ 

the same process for < t < nT\,n = 3,4, ... and finally get a global solution {g,y) £ C{B2i) 
satisfying 



\\(.0,^)\\lo^^bi,) 

+Ai| _ + y/^Vg ^ ^ \ 



< Ci||(£>o,vo)|b|, <|. (4.25) 

The choice of Si is sufficient to ensure Q+'4^ > 0. Then it follows form Remark l3.1l that (p, v) € 
C^([0,oo)xM^) is a classical solution of (|l.ip - (|1.2p with p > 0. Furthermore, we arrive at Theorem 
Owith(5o = min(Ji/2,(5i/2Ci). □ 
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5 Relaxation limit 

In this section, we give the proof of Theorem 11.21 

Proof. From (jl.7p and Remark |2. 2 1 we deduce that quantities sup^>o ||/0^ — p||_Bf ^ and 

are bounded uniformly with respect to r. Therefore, the left-hand side of (jl.4p reads as x the 
time derivative of a quantity which is bounded in L^(R"'" x R*^), plus T^X the space derivative 
of a quantity which is bounded in x W^). So, this allows us to pass to the limit r — )• in 

the sense of distributions, and we arrive at 

_ VP(/>^) ^ in P'(M+xM'^). 

T 

Inserting the weak convergence property into the first equation of (|1.4p . we have 

dsp'' - AP{p'') ^0 in V'{R+ x M*^) 

as T — 0. 

On the other hand, by (jl.4p . we conclude that dsP^ is bounded in L'^{R+ , B^i'^). Hence, it 
follows from Proposition 12.31 and Aubin-Lions compactness lemma in |18j that there exists some 
function M G C(M+, p + ^2,1) such that as r 0, it holds that 

{p^}^Af strongly in C{[0,T], {B^-^){B,)), 

for any T > and 5 £ (0, 1), which implies that A/" is a global weak solution to the porous 
medium equation ()1.6p satisfying (jl.Sp . For more details, the reader is referred to e.g. 

Therefore, the proof of Theorem 11.21 is complete. □ 



6 Appendix 

As we known, Vishik, Bahouri, Chemin and Danchin et al. j2m [H [S] have obtained some 
estimates of commutator, however, their results are unable to be applied to our case directly. 
Hence, following from their arguments, we develop a new estimate of commutator. 

Proposition 6.1. Let s>0, l<r<oo and p,pi,p2 G [l,oo]'^ with l/p = l/pi + l/p2- There 
exists a generic constant C > depending only on p,pi,p2,cr,r,d such that 

2'>^\\[f,AM9\\Lp < Cc,||V/||^,,n^.^^in^o^^^(||<7||B^.^,, + IIV^IIlpi), (6.1) 

where the operator A := div or V. As a direct consequence, when 1 < p < P2 ^ Pi ^ 00, if 

s > 1 + d( ) or s = 1 + d( ) and r = 1, 

\p2 pi/ \P2 Pi'' 

then 

2'^^\\[f,A,]Ag\\LP < Cc,\\Vf\\^s-^\\g\\Bs^^^, (6.2) 
where {cq} denotes a sequence such that \\{cq)\\ir < 1. 
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Proof. To show that the gradient part of / is involved in the estimate, we need to spht / into 
low and high frequencies: / = A_i/ + /. Obviously, there exists a constant C > such that 

||A_iV/||lp <C||V/||iP, \\Vf\\Lp<C\\Vf\\L., pG[l,oo]. (6.3) 

Since g is spectrally supported away from the origin, that is, there exists a radius < < | 
such that Supp Tff]B{0,R) = 0, Lemma [2?T] implies 

IIAgV/lliP «25||A,/||lp, pG[l,oo], q>-l. (6.4) 

Without loss of generality, we proceed the proof with Ag = divg. Taking advantage of 
Bony's decomposition, we have 

[/,Aq]div5 = [/,A,]div5+ [A_i/,Ag]div(7 

= fAgdivg - A,(/div5) + [A_i/, Ag]dwg 
= TjAgdivg + T^^divg/ + ^(/' ^9div<7) 

-A,(T^-div5 + T^^^J + R{f, divg)) + [A_i/, AJdivg. 

Set [/, Ag]dwg = ^Li where 

= TjAgdjg^ -AgTjd.g^, {divg := d,g^) 

= d,Rif,A,g^)-d,A,Rif,g^), 
Fl = A,R{djf,g^)-R{djf,A,g^) 
F^ = [A_i/,A,]div5. 

By Proposition l2.lt we have 

^5 = Y.^^'-^f^i'^'^^^a' -A,Y,Sq'-ifA,>d,g^ 

q' g' 

= Yl [Sq'-if,A,]djA,.g^ 

\q-q'\<i 

= V / h{y)[Sg,-J{x)-Sg,.J{x-2"^y)]djA,,g^ix-2-'^y)dy. 

Then, applying first order Taylor's formula. Young's inequality. Lemma |2. II and ()6.3p . we get 

2^'^IIF^^IIlp < C Y1 l|V/||LPi2('^-i)(''-^')2«''^||A,,g^||iP2 

|g-(j'|<4 

2'^'''\\Ag,g\\LP2 



< CCql\\Vf\\LPl\\g\\B-^^,, Cqi:= Yj 



q-q'\<4: "^"-°P2,r 
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and 



LP 



q'>q-3 

q'>q-3 

< C Y 2^-'^'||V/||i.i2^-||A,5||LP2 

q'>q-3 

29-||A,5||lp2 



< CCg2\\Vf\\LPl\\g\\BS^^, 0,2-.-- 



The third part is proceeded as follows: 

P'q = -\TQqlf 



then 



\q-q'\<4 



|g-q'|<4 

< C Y 2(«-'?')-||5,,_ia,9^||iPi2^'(-~i)||A,,V/||z.P2 

\q-q'\<i 

<- r llY7fll wxi II 2'?'(-i)||Ag.V/||^P. 

< Cc,3||V/||5. i||V5||lpi, c,3:= 2^ Qiivfll . . 



, , 9||V/|La 



By the definition 12.21 and Proposition I2.H we have 



= djRif,A,g^)-djA,Rif,g^) 

= Yl d,{A^,fA,,A,g^)-d,A,R{f,g^) 

k-g'l<i 

^ q ^ ^ q • 

For the first term, using (|6.4p and Lemma |2. 11 we obtain 

2''^||F^^'i||z.P < 2^^^ Y l|A,'V/||LPi||A,,5^||iP. +2'?- Y 2^-^'||A,,V/||lpi ||A,,5^- 

k-g'l<i l9-'/'l<i 

< C||V/|Upi Y 2('^-''>25'-||A,,5^||iP2 

k-g'l<i 

+C7||V/||lp, 2('?-^')(-+i)2'^'-||A,,5^|Up2 
l'2-g'l<i 

< C'Cg4(i)||V/||LPi||5||Bj^^,,, C4(i) := 2^ II • 

\q-q'\<l liyil'^P2,»- 

18 



The second term is estimated as: 

2<?-||9,A,i?(/,9^)||iP 



< 




< 


CCg4(2) 


< 


CCg4(2) 


< 


CCq4(2) 


< 


CCq4(2) 


< 


CCg4(2) 



1 1 1 

I Ri wgWB" I - = — I — 

|ri lls'lli?'^ 



|V/b0 JI^IIb^ Cg4(2) 



'Br, 



where we have used Lemma [2.21 and the result of continuity for the remainder (Proposition I2.2p . 
Among them, s + 1 > is required. 

For , it follows from the same argument as that 



2^ni^g'lUp < Cc,5||V/||iPinB0^,Jl5l|B,^^,,, 



where 



' <1 11.3 ll-Upj.,. 



4||ii(a/,g^-)||B„^ 



and s > is required. 

For F6 = ^|g_g,|<,[Ag(A_i/5jA,,5-')-A_i/A,A,,a,-5^] (g^' = by applying first 

order Taylor's formula, Young's inequality, Lemma |2. II and ()6.3p . we have 



= V / %) A_i/(x) + A_i/(x-2-'^y) 

< C ^ 2('?~«')(--i)||VA_i/|Upi2'?'n|A,,g|UP2 

k-q'l<i 

< C'Cg6||V/||LPi||5||B- , 



A,,djg^{x-2~'iy)dy 



LP 



q-q'\<l ll^ll-°P2. 



Adding above these inequalities together, the inequality (16. ip is followed with Cg = g Yli=i Ci 
satisfying ||(cg)||r < 1. 
Furthermore, if 

s > 1 + d( \ or s = 1 + ^ and r = 1 

Vp2 pi/ Vp2 Pi^ 

with 1 < p < P2 < Pi ^ OO) we have the following embedding properties: 

the inequality ()6.2p follows immediately. 

Therefore, the proof of Proposition 16.11 is complete. 



□ 
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Having Proposition 16.11 we may deal with some estimates of commutator of special form in 
the proof of a priori estimate, which are not covered by Lemma 12.31 For clarity, we give them 
by a corollary. 

Corollary 6.1. Let a = 1 + d/2. There exists a generic constant C > depending only on a,d 
such that 

( 2^-||[£,,Ag]divv||^2./.+2 <Cc,||V^||5.-i||v||B. j 

I 2^-||[v,A,]-V^?||^2./.+2 <Ccg||Vv||^.-i||i5||B.^; ^ ' > 

where {cq} denotes a sequence such that \\{cq)\\ii < 1. 

Proof. In Proposition 16.11 it suffices to take 

2d , . d , 

P = J^y^ > 2), p2 = 2, 0- = 1 + - and r = 1, 

the conclusions follow obviously. □ 

According to Holder inequality and Remark 12.21 it is ^ot difficult to achieve the estimates 
of commutators in Uj.^L?'^/'^'^'^) spaces. 

Corollary 6.2. Let a = 1 + d/2 and 1 < 9 < oo. Then there exists a generic constant C > 
depending only on a, d such that 

( 29-||[^, A,]divv||^e(^2./.+2) < ^^c,||V^>||j.i(^._i)||v||~«,(^^„^^^; 
\ 2'?-||[v,A,] • V£>||ie^(i2d/d+2) < C^Cg||Vv||~e,^^^._i^||£>||~e,^^^„^^^; 

where {cq} denotes a sequence such that ||(cq)||;i < 1 and ^ = ^ + ^• 

Remark 6.1. Actually, if we take p = p2, Pi = oo, a = 1 + d/p and r = 1 in Proposition 16.11 
we can also deduce Lemma ESI 
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